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1. Introduction: Heavy quarks on the lattice 

The violation of CP symmetry is one of the most important but still least understood 
features of the Standard Model. Much effort has been invested in order to analyse CP 
violation within the framework of the Cabibbo-Kobayashi-Maskawa (CKM) pattern of 
quark mixing. Heavy quark systems play an important role in the study of the CKM 
mixing matrix: they contain information on its least known elements and serve to test its 
parametrisation as a unitary 3x3 matrix in terms of three angles and one CP violating 
complex phase. 

The phenomcnological extraction of CKM matrix elements involving heavy quarks 
has been hampered by large hadronic uncertainties in the evaluation of current matrix 
elements appearing in the relevant weak decay amplitudes. Since quarks are confined 
within hadrons the exchange of soft gluons makes a perturbative analysis of weak decays 
impossible. During the past decade one has therefore sought to compute these matrix 
elements non-perturbatively in lattice QCD. 

Lattice simulations of QCD are by now a mature field. They allow for the evaluation 
of hadronic masses, decay constants and form factors from first principles. However, be- 
fore lattice results can be applied phenomenologically a critical assessment of systematic 
errors is required. In this lecture I shall present lattice results for leptonic and semi- 
lcptonic decays of heavy quark systems as well as B — B mixing. Systematic effects and 
methods how to increase the accuracy of lattice data will be discussed. More detailed 
information can be found in recent review articles [1—4]. 

In order to formulate QCD on a discrete grid of points one approximates space-time 
by a euclidean, hypercubic lattice with lattice spacing a and volume L 3 ■ T. One then 
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chooses a discretisation of the QCD action involving the quark fields q(x),q(x) and gauge 
fields Uft(x) <G SU(3), n = 1, . . . , 4. One such discretisation was formulated by Wilson [5] 
("Wilson fermions") and has been used for almost all results I shall describe. Using the 
discrctiscd QCD action 5qcd, one can define a partition function 



where in the last line we have integrated out the quark fields, resulting in the determinant 
of the Wilson-Dirac operator times the exponentiated pure gauge action. The expectation 
value of an observable O is defined as 



In a numerical simulation one evaluates (O) by generating a representative sample of 
N c gauge configurations in a Monte Carlo procedure. The expectation value (O) is then 
approximated by the sample average O 



where Oi is the value of the observable computed on the ith configuration. Since one 
is working with a finite number of configurations, O is obtained with a statistical error 
proportional to 1/ ^/W c . This procedure yields the observable for non-zero values of the 
lattice spacing a. The continuum result is obtained in the limit a — > 0. In practice, this 
usually means that the simulation must be repeated for several values of a so that the 
results can be extrapolated to a = 0. 

The appearance of the determinant in eq. (1) presents a major obstacle for numerical 
simulations since it is highly non-local and thus its evaluation is a huge computational 
overhead. Before the advent of efficient algorithms it had been suggested to disregard 
its effects completely by setting det Q = 1. This defines the so-called Quenched Ap- 
proximation [6], which in physical terms corresponds to neglecting quark loops in the 
determination of (O) . Although of unknown quality, the quenched approximation is still 
widely used and also accounts for most of the material covered here. 

Since the (inverse) lattice spacing acts as an ultraviolet cutoff, its value in physical 
units places constraints on the scales that one is able to study. On present computers 
typical values lie in the range of a -1 = 2 — 4GeV. These relatively low values of the 
cutoff imply that one expects large cutoff effects for the charm quark whose mass is not 
too far below a -1 [GeV]. More importantly, the b quark cannot be studied directly since 
its mass lies above the cutoff. The following methods are used used to circumvent this 
problem: 

• Reduction of lattice artefacts 



(1) 



Z = 




(2) 




(3) 
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• Static approximation 

• Non-relativistic QCD (NRQCD) 

In the first approach one seeks to cancel the leading cutoff effects of order a in the Wilson 
action by employing so-called improved actions and operators [7—11]. For instance, the 
0(a) improved lattice action is defined as 

(4) SQ CD [U,q,q] = S Q cn[U,q,q] +c sw ^ ^ q{x)a llv F llv {x)q{x), 

where c sw is an improvement coefficient and F^ u is a lattice transcription of the field 
tensor. Further coefficients appear in the definitions of the improved axial and vector 
currents. Provided that all improvement coefficients are chosen appropriately one can 
show that lattice artefacts of 0(a) are cancelled completely in masses and matrix el- 
ements. In a different approach, which can be applied for improved and unimproved 
actions, a reduction of lattice artefacts is achieved through absorbing higher-order effects 
of the quark mass into a rescaling factor [12] 

(5) q(x) — ► e amp/2 q(x), am P = ln(l + am,), 

where am q is the bare quark mass. This normalisation is called the El-Khadra-Kronfeld- 
Mackenzie (EKM) norm. It must be emphasised that neither improvement nor the EKM 
norm solve the problem that the b quark cannot be studied directly. However, they 
reduce lattice artefacts around the charm quark mass, so that the obtained results can 
be extrapolated to the b quark mass much more reliably. 

The static approximation is based on the leading term of the expansion of the heavy 
quark propagator in powers of the inverse heavy quark mass 1/mg. Thus, one regards 
the b quark as infinitely heavy in this approach and expects that results will be subject 
to corrections of order Aqcd/^q- Finally, NRQCD is an effective theory based on an 
expansion of the QCD action in the four-velocity of the heavy (non-relativistic) quark. 
Again one expects corrections to the effective theory whose influence on physics results 
has to be assessed. 

From this discussion it is clear that none of the above methods gives an entirely 
satisfactory description of heavy quarks on the lattice. However, they all provide com- 
plementary information which can be used to reveal the full picture. 

Besides lattice artefacts, another important source of systematic errors is the explicit 
breaking of chiral symmetry by the Wilson action. As a consequence, lattice versions of 
the local vector and axial currents are not conserved. Instead they are related to their 
continuum counterparts by finite renormalisations Zy and Za, respectively. Although 
Zy and Za have been computed non-perturbatively for 0(a) improved actions [13,14] 
these factors are known only in one-loop perturbation theory in the unimproved case. 
Furthermore, explicit chiral symmetry breaking causes operators with definite chirality - 
such as the four-fermion operator used to describe B — B mixing - to mix with operators 
of opposite chirality. Therefore, several matrix elements have usually to be determined 
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on the lattice and subsequently matched to the continuum matrix element. The general 
expression in the unimproved theory thus is 

(6) (f\0\ir^ = J2z a {f\0^\i) + 0(a), 

a 

where the Z a 's are the appropriate normalisation factors, and lattice artefacts of order a 
arise through mixing with higher dimension operators. 

Finally lattice estimates of dimensionful quantities are subject to uncertainties in 
the lattice scale. They are due to the fact that different quantities like f v , M p , . . ., which 
are commonly used to set the scale a -1 in physical units give different results. This is 
closely related to using the quenched approximation, since loop effects are not expected 
to be the same for different quantities. 



2. Leptonic decays of heavy-light mesons 

The leptonic decay constant fp of a heavy-light pseudoscalar meson is related to the 
matrix elements of the axial current on the lattice via 



(7) <0|Ai(0)|PS) = Mpf P /Z A , 

where Mp is the pseudoscalar mass and Za is the renormalisation factor of the lat- 
tice axial current. Both the matrix element and Mp are obtained from the asymptotic 
behaviour of the euclidcan correlation function of the axial current at large separation t 

(8) 5> 4 (£, i)4(0)> ^° l(0|A4 o ( ° )|PS)|2 f e-*** + e-^(T-t)) . 

— ' 2M p L J 

X 

The decay constant fp can then be studied at several different values of the mass of the 
heavy quark. This enables one to study some predictions by the Heavy Quark Effective 
Theory (HQET). It is well known that HQET predicts the following scaling law in the 
limit of an infinitely heavy quark 



(9) fpVMp ^ const x a s (M P ) 



-VPo 



where a s is the strong coupling constant and (5q = 11 — 2n//3. In order to test the 
quality of this prediction we plot in Fig. 1 the quantity 



(10) *(Mp) = f P 




a s (M P ) 
a s (M B ) 



2//3o 



as a function of 1/Mp, using data in the static approximation [15-17] and for relativistic 
heavy quarks [16,18,19]. Using the static approximation as the limiting case, the figure 
illustrates that there are large corrections in 1/Mp to the scaling law, provided that 
lattice artefacts have been treated, cither by using improvement (i.e. c sw > 1) or by 
employing the EKM norm. Failure to address the problem of lattice artefacts leads to an 
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Fig. 1. - $(Mp) versus 1/Mp using data from different collaborations for a ~ 3 GeV. The 
data from [19] are shown with (plus signs) and without (crosses) including the EKM factors. 



inconsistent mass behaviour of the data in the static and relativistic regimes, (c.f. crosses 
in Fig. 1). We conclude that the HQET scaling law is not satisfied at the physical B and 
D meson masses and that the treatment of lattice effects is crucial for the computation of 
heavy-light decay constants in general. Moreover, this example illustrates the interplay 
between different formulations of heavy quarks. 

Recent results for heavy-light decay constants using relativistic heavy quarks and 
NRQCD are compiled in Table I. Besides the statistical error, most groups now quote 
one or more systematic errors. Note, however, that the estimation of systematic errors 
can vary significantly among different collaborations. The results shown in Table I are 
broadly consistent, but the errors are still large and dominated by systematic effects. 
Despite the apparent consistency of the data it would be premature to simply combine 
them into global estimates of decay constants, because the details of the analysis of the 
raw lattice data can differ substantially among different groups. A common analysis of 
lattice data from many collaborations has been described in [4]. There the aim was to 
perform the extrapolation to the continuum limit and to present a uniform estimation of 
systematic errors from a number of sources. These include the choice of lattice scale (e.g. 
fn, M p , . . .), the quark field normalisation, uncertainties in the perturbative values of 
Z\, and variations in fitting and extrapolation procedures. The results in the continuum 
limit can be summarised as 



(11) f D = 191 ± 19 (stat) t 2 l (syst) MeV, f Ds = 206 ± 17 (stat) + 2 | (syst) MeV 
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Collab. 



a [fm] f D [MeV] 



f Ds [MeV] 



f B [MeV] 



f Bs [MeV] 



FNAL* [20] 





205(9)(27) 


215(7)(30) 


166(10)(28) 


1.17(4)(3) 


JLQCD* [21] 





192(10) 11£ 


213(11) ±» 


163(12) til 




MILC* [22] 





186(10) ±U±° 


199(8) it? - 1 " 


153(10) tT 3 t 13 


1.10(2) ±1 1% 


APE [23] 


0.07 


221(17) 


237(16) 


180(32) 


1.14(8) 


LANL [24] 


0.09 


229(7) t? 6 


260(4) til 






PCW [19] 





170(30) 




180(50) 


1.09(2)(5) 


UKQCD [18] 


0.07 


185+4 t 4 7 2 


212(4) t 4 r 


160(6) til 


1.22 tt 


BLS [16] 


0.06 


208(9) (35) (12) 


230(5)(10)(19) 


187(10) (34) (15) 


1.11(6) 


Hirosh. 1 [25] 


0.12 






184(7)(5)(37)(37) 


1.23(3)(3) 


GLOK 1 " [26] 


0.09 






183(32)(28)(16) 


1.17(7) 


SCO 1 [27] 


0.10 






126-166 


1.24(4)(4) 



Table I. - Results for heavy-light decay constants from different collaborations. Data marked by 
an asterisk are preliminary. Results obtained using NRQCD are marked by a dagger. All other 
collaborations have used relativistic heavy quarks. The convention f^ — 131 MeV is understood. 



(12) f B = 172 ± 24 (stat) t\l (syst) MeV, f B Jf B = 1-14(8) 

Further details and comparisons to other theoretical results can be found in [4] . 

3. — B — B mixing 

We now turn our attention to matrix elements which arc used to describe oscillations 
between B° and B states. The mechanism of B° — B mixing is illustrated by the box 
diagrams shown in Fig. 2. 

, u, c, t 

o — ► 




u, c, t 

Fig. 2. - Box diagram contributions to B° — B° mixing. 




The mass difference Am between the mass eigenstates B° and B is related to the 
CKM matrix elements V t d and V t b via 

(13) Amcx \V td V; b \ 2 M B f B B B . 

Here, f B is the decay constant of the B meson encountered in the previous section, and 
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Bb denotes the B parameter defined by 



(14) 



B B (li) = 



(B°\O L (ri\B ) 



3-1 B 1V± B 

where \i is the renormalisation scale and the four-fermion operator Ol is given by 
(15) O l = (67„(1 - 76)d) (&7^(1 - 76)d) . 

The "hat" on £?b in eq. (13) signifies that the dependence of Bb on the renormalisation 
scale has been divided out. The resulting renormalisation group invariant B parameter 
can be defined at leading (LO) or next-to-leading order (NLO) in a s via 



(16) 
(17) 



B B ° = a s (»)-W°B B (ri 



B^° = a s (^) 



-2/0o 



1 



47T 



where J5 is derived from the one- and two-loop anomalous dimensions of the operator 
Ol- It is important to realise that the combination J b Bb is the principal unknown 
quantity which relates the experimentally measured mass difference Am to the CKM 
matrix elements in eq. (13). Thus, f B B B is an important ingredient for the study of CP 
violation. 




* B B {n) 





7 2 
Z L 



Fig. 3. - Lattice measurable for the contribution of Ol to the B parameter (see text). Full dots 
represent insertions of the axial current. 



On the lattice Bb(h) is obtained from a ratio of two- and three-point functions 
as depicted in Fig. 3. Apart from the renormalisation constant Za of the lattice axial 
current the relation between the lattice and continuum versions of Bb will also contain 
the factor Z L associated with the lattice version of Ol- Furthermore, as there is mixing 
between Ol and other four-fermion operators with opposite and mixed chiralitics, all 
contributions have to be evaluated and matched to the continuum theory. Since the rel- 
evant Z-factors have only been computed in one-loop perturbation theory, the matching 
procedure introduces considerable uncertainties into the final estimates of the B param- 
eter. For lattice results obtained in the static approximation and using c sw = 1, these 
uncertainties have been estimated to be as large as 25% [4, 17,28]. 
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Collab. 


a [fm] 


B B {m b ) 


nLO 




Kentucky 1 " [30] 


0.09 


0.97(4) 


1.42(6) 


1.54(6) 


G+M 1 " [28] 


0.09 


0.63(4) 


0.92(6) 


1.00(6) 




0.09 


0.73(4) 


1.07(6) 


1.16(6) 


UKQCD 1 " [17] 


0.07 


0.69 tl -i 


1.021* t% 


1-10 ±1 ±1 


UKQCDt [29] 


0.07 


0.83 tl tl 




1.32 1* tl 


JLQCD* [31] 


0.06 


0.840(60) 


1.23(9) 


1.34(10) 




0.08 


0.895(47) 


1.31(7) 


1.42(8) 


B+S* [32] 





0.89(6) (4) 


1.30(9)(6) 


1.42(10)(6) 


ELC [33] 


0.05 


0.84(5) 


1.24(7) 


1.34(8) 


BDHS [34] 


0.08 


0.93(14) 


1.36(20) 


1.48(22) 



Table II. - Data for the B parameter from different collaborations at a reference scale m b = 
5GeV. Data marked by a dagger have been obtained using the static approximation. Data 
marked by an asterisk are preliminary. 



Lattice data for Bb(p) have been published for propagating and static heavy quarks 
and are shown in Table II. It is clear that lattice data for Bb do not yet allow for a 
continuum extrapolation as in the case of fs- Instead one may quote a common estimate 
with an error that encompasses the spread of different results. In ref. [4] the result is 

(18) B B (5 GcV) = 0.85 til, B B = 1.3 tj- 
The global result for the ratio BbJBb quoted in [4] is 

(19) B B JB B = 1.00 ±0.02. 

We can now combine the results for the decay constant /b in eq. (12) with that for the 
B parameter. Combining the errors in quadrature one finds 

(20) Ib\[b~b = 195ll°MeV. 

For the SU(3)-flavour breaking ratio involving both the decay constants and B parameters 
the global result in [4] is 



fl 3 B B B 
S%Bb 

These estimates can now be used in the study of the CKM matrix and CP violation. 



(21) = i.3 8 ±0.15. 

Jb b b 



4. Semi-leptonic B — > it and B — > p decays 

We will now discuss lattice results for semi-leptonic decays of B mesons. Due to 
lack of space we will concentrate on heavy-to-light transitions. Semi-leptonic B — > D 
and B — > D* decays are reviewed, for instance, in ref. [1]. 
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Decays like B° — > tt + 1~vz or B° — > p + l~vi have attracted much interest recently, 
since they can be used to extract V u b, which is one of the most poorly known CKM 
matrix elements. Compared to the pseudoscalar decay constant or the B parameter, the 
study of scmi-lcptonic decays is more complicated due to the kinematics involved in their 
description. The relevant matrix elements of the weak V — A current are parametrised 
in terms of form factors which depend on the momentum transfer q 2 between the initial 
and final mesons. For instance, if the final state is a pseudoscalar meson only the vector 
current contributes, and there are two independent form factors /+ and f° 

(22) (PS(fc) K\ B(p)) = /+(«») {(P + fc)„ - + /V)^^, 

where = — fc M . If the final state is a vector meson, there are four form factors 
V(q 2 ), A 1 {q 2 ) 1 A 2 (q 2 ) and A (q 2 ), associated with matrix elements of the vector and 
axial vector currents, respectively. 

The need to control lattice artefacts places constraints on the possible values of lat- 
tice momenta p and k. Together with the constraints on the heavy quark mass this 
implies that one usually obtains form factors for typical momenta \p\ < 1.5GeV/c 
and heavy quark masses in the region of charm. The "generic" heavy-to-light semi- 
lcptonic decay in lattice simulations is thus D — > Kdvi for momentum transfers in the 
range — 1 GeV 2 /c 2 < q 2 < 2GeV 2 /c 2 . Although in principle one is interested in the 
<7 2 -dependence of form factors, they are commonly quoted at q 2 = 0. Frequently a simple 
pole ansatz is used to model the g 2 -dependence: 

(23) F(q 2 ) = F{0) 



(l-9 2 /M 2 olc )"-' 

Here, F is a generic form factor, M po i c is of the order of the heavy-light meson mass, and 
riF — 0, 1, 2, . . . parametrises constant, monopole, dipole and higher multipole behaviour 
of F. A recent compilation of lattice results for form factors for semi-leptonic D decays 
can be found in [1,35]. Since the lattice form factors for these decays are determined in a 
region around q 2 = 0, the pole ansatz in eq. (23) merely serves to guide the interpolation 
of F(q 2 ) to q 2 = 0, so that no model dependence is introduced. 

The situation changes significantly if one considers semi-leptonic B — > 7r or B — > p 
decays. Here the form factors are obtained through extrapolation in the heavy quark mass 
to the mass of the b quark. Since similarly large values of lattice momentum \p\ cannot be 
considered due to restrictions imposed by lattice artefacts, the accessible region of q 2 is 
pushed to large values near q^^. This in turn leaves a long and potentially uncontrollable 
extrapolation in q 2 to determine F(0). In order to map out the g 2 -behaviour one usually 
cannot avoid relying on model assumptions. 

We shall now describe how the <7 2 -behaviour can be constrained by requiring con- 
sistency with heavy quark symmtery (HQS), kinematical constraints and scaling laws 
implied by light-cone sum rules. In analogy to the decay constant in eq. (9), HQS pre- 
dicts the following leading scaling behaviour of form factors in the infinite mass limit at 
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Collab. 


a [fm] 


/ + (0) 


V(Q) 


MO) 


A 2 (0) 


UKQCD [38] 


0.07 


0.27(11) 


0.35 1!? 


0.27 ±1 


0.25 ±1 


GSS [39] 


0.06 


0.43(19) 


0.65(15) 


0.28(3) 


0.46(23) 


APE [40] 


0.09 


0.35(8) 


0.53(31) 


0.24(12) 


0.27(80) 


ELC [41] 


0.05 


0.30(H)(5) 


0.37(11) 


0.22(5) 


0.49(21)(5) 



Table III. - Lattice results for form factors for semi-leptonic B — > tt and B — + p decays. 



fixed values of w (which is the product of four- velocities of the initial and final mesons) : 

(24) /+H-M 1 / 2 , /°(w)-M- 1 / 2 , V(uj)^M^ 2 , Ai(u) <~ M~ 1/2 , . . . 

Here, M is the mass of the heavy-light meson. The extrapolations of form factors to the 
b quark mass can be performed using model functions motivated by the above scaling 
laws. Additional scaling laws are provided by light-cone sum rule analyses [36,37], which 
predict that all form factors scale like F ~ M~ 3 / 2 at q 2 = 0. In the heavy quark limit 

(25) '-sfc^ 

and thus the scaling laws predicted by HQS and light-cone sum rules can be combined to 
infer a value for rip in the pole formula eq. (23). Also, kinematical constraints at q 2 = 
such as 

(26) /+(0) = /» 

can be used to analyse the ^-dependence of lattice form factors. In Table III we list 
lattice estimates for form factors for B — > n + t~Vi> and B — ► p+Crv^ decays. It should 
be noted that only the UKQCD data are consistent with all constraints discussed above. 

A very different approach to constrain the g 2 -dependence in a modcl-indcpcndcnt 
fashion has been discussed by Lcllouch [42]. Lattice data for the form factors / + and 
/° for B — > ir + l~Vit obtained near have been combined with dispersion relations 
and the kinematical constraint eq. (26). The method relies on perturbative QCD in the 
evaluation of the dispersion relations and general properties such as unitarity, analyticity 
and crossing. However, existing lattice data for the form factors are at present not precise 
enough in order to allow for stringent bounds at q 2 = 0. 

Another proposal to avoid model dependence in the extraction of V u t, was made 
in [43]. Here one concentrates on the exclusive decay B — > p + £~17e in the region near 



q^x- Instead of attempting to extract the form factors at q 2 = one parametrises the 
differential decay rate by 



(27) ^ = io- 12 §^ 2 VMF)A\i + B(q 2 - ,LJ), 
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where A and B are parameters and A is a phase-space factor. The combination A 2 {I + 
B (q 2 — g^ax)) parametrises the long-distance hadronic dynamics, and A 2 provides the 
overall normalisation. Using lattice data for the form factors to evaluate the differential 
decay rate the authors of [43] find 

(28) ^ = 4.6lg;|±0.6GeV, B = (-8 ±\) • 1(T 2 GeV 2 . 

The corresponding prediction of the decay rate is shown in Fig. 4. Given sufficient ex- 
perimental data for the decay rate in conjunction with accurate lattice results, a deter- 
mination of V u b will be possible. 

* * * 

I am grateful to Profs. I. Bigi and L. Moroni for the kind hospitality at this school. I 
would like to thank J. Flynn, L. Lellouch, J. Nieves, C. Sachrajda and other members of 
the UKQCD Collaboration. The support of PPARC through the award of an Advanced 
Fellowship is gratefully acknowledged. 
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